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Abstract. The paper deals with a 3-parameter family of probability measures on 
the set of partitions, called the z-measures. The z— measures first emerged in connec- 
tion with the problem of harmonic analysis on the infinite symmetric group. They 
are a special and distinguished case of Okounkov's Schur measures. It is known that 
any Schur measure determines a determinantal point process on the 1— dimensional 
lattice. In the particular case of z— measures, the correlation kernel of this process, 
called the discrete hypergeometric kernel, has especially nice properties. The aim of 
the paper is to derive the discrete hypergeometric kernel by a new method, based 
on a relationship between the z— measures and the Meixner orthogonal polynomial 
ensemble. In another paper (Prob. Theory Rel. Fields 135 (2006), 84-152) we apply 
the same approach to a dynamical model related to the z-measures. 



Introduction 

Main definitions and motivations. Recall that a partition is an infinite mono- 
tone sequence of nonnegative integers, A = (Ai > A2 > • ■ ■ ), with finitely many 
nonzero terms Ai . There is a natural identification of partitions with Young dia- 
grams; for this reason, we denote the set of all partitions by symbol Y. Clearly, 
Y is a countable set. To each partition A 6 Y we assign a weight depending on 
three parameters z, z', and £. Under suitable restrictions on the parameters (for 
instance, if z and z' are complex numbers conjugate to each other and < £ < 1) 
all the weights are nonnegative and their sum equals 1. Then we get a probability 
measure on the set Y, which makes it possible to speak about random partitions. 
The measures on Y obtained in this way are called the z-measures and denoted as 
M z>z i£ (see section 1 for precise definitions). 

Our interest in the z-measures is mainly motivated by the fact that they play 
a crucial role in harmonic analysis on the infinite symmetric group, see [KOV1], 
[KOV2], [B02], [01]. On the other hand, for special values of parameters z,z' the 
z-measures turn into discrete orthogonal polynomial ensembles which in turn are 
related to interesting probabilistic models: the directed percolation model [Jol], 
the stochastic growth model of [GTW] , random standard tableaux of rectangular 
shape [PR]. The z-measures are studied in many research papers: [B02], [B03], 
[B04], [B05], [B06], [BOS], [Ok2] ; see also the expository papers [BOl], [01]. 
Finally, note that the z-measures are a particular case of more general objects, the 
Schur measures introduced by Okounkov in [Okl] and further investigated by many 
people. 

Although the z-measures are quite interesting by themselves, the main problems 
concern their limits as parameter £ approaches the critical value 1 (parameters z, z' 
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being fixed). Note that, as £ — > 1, the weight of each partition tends to 0, that 
is, the measure runs away to infinity. Thus, to catch possible limits we have to 
embed Y in a larger space. It turns out that there are different limit regimes, 
and for each regime the limit measure lives on a suitable space of infinite point 
configurations (see our paper [B05] for more details). In other words, the limit 
measure determines a random point process. An appropriate way to describe point 
processes is to use the language of correlation functions, and the first necessary step 
is to interpret the initial z-measures as point processes, too. 

To this end, we use a well-known interpretation of partitions as Maya diagrams, 
which are semi-infinite point configurations on the 1-dimensional lattice. It is 
convenient to identify the lattice with the subset T! := Z + \ C M of (proper) 
half-integers. Then the Maya diagram of a partition A e Y is the configuration (or 
simply the subset) {A^ — i + | | i = 1, 2, . . . } C U . Each z-measure M zz >£ thus 
gives rise to a random point configuration on Z' (or a point process on Z'), and its 
nth correlation function p n (n = 1, 2, . . . ) expresses the probability p n (x\, . . . , x n ) 
that the random configuration contains an arbitrary prescribed finite set of points 
Xi, . . . ,x n in Z' . 

It is worth noting that the correlation functions survive in various limit regimes, 
which explains their efficiency. 

A remarkable property of the z-measures is that, for any n — 1,2,..., the prob- 
ability p n {x\, . . . , x n ) can be written as the n x n determinant det[K(x i} xj)] where 
K{x, y) is a function on U x Z' not depending on n (it depends on parameters z, z' , £ 
only). Random point processes with such a property are called determinantal ^ and 
the function K(x, y) is called the correlation kernel. 

As was first shown in [B02], the correlation kernel of the z-measure M ZjZ i £ can 
be explicitly written in terms of the Gauss hypergeometric functions; for this reason 
we called it the discrete hypergeometric kernel. Then a number of different proofs 
were suggested in [Okl] (see also [BOk]), [Ok2], [BOS]. The goal of the present 
paper is to better understand the nature of this kernel. 

The results. Now we are in a position to describe our main results: 

(1) We introduce a system of functions ip a (x) = ipaix', z, z', £), where the triple 
(z, z 1 , £) is the parameter of the z-measure, x is the argument ranging over Z', and 
a e U is an additional parameter. For fixed (z, z' , £) and varying a, the family {ip a } 
forms an orthogonal basis in the coordinate Hilbert space £ 2 (Z'). Each function tp a 
can be expressed through the Gauss hypergeometric function. 

(2) We exhibit a second order difference operator D = D(z,z',£,) on U which 
is diagonalizcd in the basis {ip a }- The eigenvalue of D corresponding to the eigen- 
function ip a is equal to a(l — £). (We assume < £ < 1, so that the eigenvalue 
a(l — £) is positive or negative depending on the sign of parameter a € Z'.) 

(3) Set Z' + = {§, §,§,■■■}■ We prove that the discrete hypergeometric kernel 
can be written as 

K(x,y) = ^a{x)i>a{y), 
ae2' + 

which means that K(x,y) is the kernel (or simply the matrix) of the spectral pro- 
jection operator in £ 2 (Z') corresponding to the positive part of the spectrum of 



1 This term, introduced in [B02] and then employed in Soshnikov's expository paper [S], is now 
widely used. Earlier works used the term "fermion point processes" . 
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D. This formula together with a three-term recurrence relation satisfied by the 
eigenfunctions ip a implies another expression for the kernel: 



Thus, K(x,y) is a discrete integrable kernel (see [B] for the definition). 

(4) The above sum expresses the kernel as a series of products of hypergeometric 
functions. On the other hand we can represent the kernel by a double contour 
integral involving elementary functions only. 

From these claims one can readily derive all known results concerning the discrete 
hypergeometric kernel. 

The method. Our approach relies on the observation made in [B02] which relates 
the z-measures to the Meixner orthogonal polynomials. Assume one of the param- 
eters (z, z') is a positive integer N — 1,2,... while the other parameter is a real 
number greater than N — 1. This is a rather special degenerate case: the weight as- 
signed to a partition A vanishes unless Xn+i = = • • • = 0, so that the relevant 
partitions A depend only on the first N coordinates Ai, . . . , Ajv- It turns out that 
in this case the random TV-point configuration {Ai + N — 1, A 2 + N — 2, ... , \ N } 
on the set Z + of nonnegative integers is a well-known object: it is an example of 
a (discrete) orthogonal polynomial ensemble. The orthogonal polynomial ensem- 
bles were extensively studied in connection with random matrix theory as well as 
various discrete probabilistic models, see [De2], [Jo2], [Jo3], [K6]. In particular, 
it is well known that they are determinantal processes and their correlation ker- 
nels are closely related to the Christoffcl-Darboux kernels for the corresponding 
family of orthogonal polynomials (in our situation these are the classical Meixner 
polynomials). 

The idea of our approach to the z-measures is to regard them as the result 
of an analytic continuation of the Meixner orthogonal polynomial ensembles in 
parameter N. In particular, our difference operator D on the lattice U comes 
from the Meixner difference operator on Z+. It is worth noting, however, that the 
procedure of analytic continuation is rather delicate, because we extrapolate from 
the discrete values z = N = 1, 2, . . . to continuous values z e C. It is this analytic 
continuation procedure that we regard as the main achievement of the present 
paper. Even though we use it to rederive a known result, in a more complicated 
dynamical situation (see the next paragraph) this method is crucial for obtaining 
new results. 

Note that instead of the Meixner polynomials one could equally well use the 
Krawtchouk orthogonal polynomials (see section 4) . 

Note also that analytic continuation of a correlation kernel off the integral values 
of a parameter was used in [Ni] in a very different situation. We are grateful to the 
referee for this remark. 

Dynamics. The present paper can be viewed as an introduction to our paper 
[B07] where the same approach is applied to studying a dynamical model related 
to the z-measures. There we derive a dynamical (i.e., time-dependent) version of 
the discrete hypergeometric kernel, K(s, x; t, y), where x and y are, as before, points 
of the lattice Z' while s and t are time variables. We also evaluate the asymptotics 
of the kernel K(s,x;t,y) in two limit regimes. We refer to [B07] for more details. 



K(x,y) 



Vii 7 ? ^-^(x)ipi(y) - ipi(x)ip_i(y) 
l-£ x-y 
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The difference operator D introduced in the present paper plays an important 
role in the dynamical picture, too. We regard this operator as the key to under- 
standing the nature of the point processes connected to the z-measures. 

Plancherel measure. In the limit as z and z' go to oo and £ goes to in such 
a way that the product converges to a positive number 9, the z-measure 
M ZtZ ',£ turns into the poissonizcd Plancherel measure Mg with Poisson parameter 
9. Our results about the correlation kernel have counterparts for Mg, see [BOO]. 
A dynamical model related to Mg is studied in [BO 8]. 

Organization of the paper. In section 1 we recall the definition of the z- 
measures and explicitly describe their relationship to the Mcixner orthogonal poly- 
nomials ensembles. In section 2 we introduce the difference operator D and we 
study in detail its eigenfunctions ip a {x). In section 3 we compute the correlation 
kernel. In section 4 we briefly discuss the relationship between the z-measures and 
the Krawtchouk orthogonal polynomial ensembles. 

Acknowledgements. Both authors were partially supported by the CRDF grant 
RUM1-2622-ST-04. The first author (A. B.) was also partially supported by the 
NSF grant DMS-0402047. 

1. Z-MEASURES 

As in Macdonald [Ma] we identify partitions and Young diagrams. By Y„ we 
denote the set of partitions of a natural number n, or equivalently, the set of Young 
diagrams with n boxes. By Y we denote the set of all Young diagrams, that is, 
the disjoint union of the finite sets Y„, where n = 0, 1, 2, . . . (by convention, Yo 
consists of a single element, the empty diagram 0). Given A € Y, let |A| denote the 
number of boxes of A (so that A G Ym), let £(X) be the number of nonzero rows in 
A (the length of the partition), and let A' denote the transposed diagram. 

By dim A we denote the number of standard tableaux of shape A. A convenient 
explicit formula for dim A is 

dimA = nj v ' A '' „ II (A<-*-A,+j), (1.1) 

where N is an arbitrary integer > £(X) (the above expression is stable in N). 
We shall need the generalized Pochhammer symbol {z)\: 
t(X) 

{z) x = Y[(z-i + l) Xi , zeC, A e Y, 



where 

(x)k = x(x + 1) . . . (x + k - 1) 



r(a; + fc) 

-Y[xT 



is the conventional Pochhammer symbol. Note that 

(*)a= n (*+j-o 

(product over the boxes of A), which implies at once the symmetry relation 

(z) A = (-l)l A l(-z) v . (1.2) 

Obviously, if z = N, where N = 1, 2, ... , then (z)\ vanishes for all A with £(A) > N. 
Likewise, if z = —N then (z)\ vanishes when £(X') = \\ > N. 
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Definition 1.1. The z-measure with parameters z, z', and £ is the (complex) 
measure M ZjZ i£ on the set Y which assigns to a diagram A e Y the weight 

M^(A) = (i-r'e |A| (m (iaf) 2 • (L3) 

The above expression makes sense for any complex z, z' and any £ G C\ [1, +oo): 
Indeed, we may assume — n < arg(l — £) < ir and then we set 

(I _ £) zz ' = |1 — £| zz ' e izz ' ar s( 1_ 4) 

Note that the weight is invariant under transposition z <-» z'. Note also the 
symmetry relation 

M_ 2 ,_^, € (A) = M^,, € (A'), 

which readily follows from (1.2). Finally, note that the z-measures are a particular 
case of the Schur measures introduced in [Okl] . 

Proposition 1.2. 1/ < £ < 1 and parameters z,z' satisfy one of the three 
conditions listed below, then the z-measure is a probability measure on Y. 

The conditions are as follows. 

• Principal series: The numbers z, z' are not real and are conjugate to each 
other. 

• Complementary series: Both z, z' are real and are contained in the same open 
interval of the form (m, m + 1), where m e Z. 

• Degenerate series: One of the numbers z, z' (say, z) is a nonzero integer while 
z' has the same sign and, moreover, \z'\ > \z\ — 1. 

Proof. As follows from [B05, §1], the series M Z ^ Z ^(X) absolutely converges and 
its sum equals 1 for any complex z,z' and any complex £ with |£| < 1. Thus, 
it suffices to check that the weights are nonnegative under the assumptions listed 
above. Since £ € (0, 1), this means that the product {z)\{z')\ is nonnegative. 

For the principal series, (z)\ and (z')\ are conjugate to each other and do not 
vanish, and for the complementary series these are both real numbers of the same 
sign. Thus, their product is always strictly positive. 

Examine now the case of the degenerate series. Assume z = N = 1, 2, . . . and 
z' > N — 1. If £(X) < N then both (z)\ and (z')\ are strictly positive, and if 
> N then (z)\ = so that the weight vanishes. Likewise, if z — —N and 
z' < —(N — 1) then the weight is strictly positive if ^(A') = Ai does exceed N, and 
vanishes otherwise. □ 

From now on we assume that the z-measure belongs to one of these three series 
and is, therefore, a probability measure. Consequently, we may speak about random 
Young diagrams, with reference to the z-measure. 

As it is seen from the above proof, for the principal series or the complementary 
series, the support of the z-measure is the whole set Y, while for the degenerate 
series, the support is a proper infinite subset of Y. 

In the remaining part of the section we will describe the relationship between the 
degenerate scries and the Meixncr polynomials. We start with a general definition. 
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Definition 1.2. Let X be a discrete subset of K, finite or countable, and let W(x) 
be a positive function on X. The TV-point orthogonal polynomial ensemble with 
weight function W is the random TV-point configuration in X such that the prob- 
ability of a particular configuration x\ > ■ ■ ■ > xm 1 where x\, . . . , xjv € X, is given 

by 

N 

Prob{xi, . . . , xn} = constjv Y\_ W(xi) ]~j (xj — x j) 2 ■ (1-4) 

i=l l<i<j<Af 

Here we assume that the cardinality of X is no less than N and that 

f " 1 
const^ 1 := sJI^^*) II {xi - Xj) 2 \ < +oo. 

xi> — >xn I i=l l<i<j<N J 

For finite X, this condition is trivial, and for infinite X, it just means that the weight 
function W has at least N — 1 finite first moments. 

The term "orthogonal polynomial ensemble" is related to the following well- 
known fact. Let Po = l,P\,Pi,... be the orthogonal polynomials with weight 
function W, the result of Gram-Schmidt orthogonalization of l,x,x 2 ,... in the 
weighted £ 2 space £ 2 (X, W). Denote by K^{x,y) the JVth Christoffel-Darboux 
kernel multiplied by y/W(x)W(y): 



N-l 



K N (x,y) = ^W(x)W(y)J2 



i=0 



Pi{x)Pi(y) 
WPW 2 ' 



where the norm refers to the weighted I 2 Hilbert space I 2 (X, W) . Note that the 
kernel Kn corresponds to the projection operator in l 2 (3C, W) whose range is the 
linear span of 1, x, ... , x^ -1 . Then we have 

Proposition 1.3. The probability that the random N -point configuration, as spec- 
ified in Definition 1.2, contains a given n-point set {yi, ■ ■ ■ ,y n } C X equals the 
determinant of the nxn matrix \K N {y i ,yj)]. 

For a proof, see, e.g., [De2], [K6, Lemma 2.8]. Note that the determinant auto- 
matically vanishes if n > N, because the kernel has rank N. 

We will be dealing with a concrete example of the weight function. This is the 
Meixner weight function, which is defined on the set X = Z + := {0, 1 , . . . }, depends 
on parameters (3 > and £ G (0, 1), and is given by 

T^Mcixncrz-N _ {P)xj X _ T(f3 + x)£ X 

W ^ {X) - x! - T{P)x\ ' G + 

(we denote a point of Z + by x instead of x because this notation is used below in 
§§2-3). 

For N = 1, 2, ... , let Y(JV) C Y denote the set of diagrams A with £(\) < N. The 
following correspondence is a bijection between diagrams A € Y(JV) and iV-point 
configurations on Z + : 

A< — ► {xi, . . . ,x N }, Xi = Xi + N-i (i = l,...,N). (1.5) 



The next fact was pointed out in [B02]: 
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Proposition 1.4. Under correspondence (1.5), the z-measure of the degenerate 
series with parameters [z = N, z' = N + [3 — 1, £), where (3 > and £ £ (0, 1), 
turns into the N-point Meixner orthogonal polynomial ensemble with parameters 
(AC)- 

Proof. It suffices to check that if £(X) < N and {xi, . . . , xn} is given by (1.5) then 
the right-hand side of (1.3) can be written as the right-hand side of (1.4) with the 
Meixner weight function. 
By virtue of (1.1), 

_ rii<j(^i _ &j) 



ia|! ) m^-r ' 

Next, with z = N and z' = N + (3 — 1 we have 



and 



N 



= £-N(N-l)/2 . 

i=l 

Combining these formulas we get (1.4) with W = W^eixner an( j 



N 



COnsW = (1 - ^N(N+P-l) r N(N-l)/2jJ 



f = \r(N -i + l)T(N -i + [3) ' 
□ 

Remark 1.5. Let A be the random Young diagram distributed according to a 
z-measure M z . z >^. Then the number of boxes |A| has the negative binomial distri- 
bution on Z + with parameters zz' and £: 

Prob{|A| = n} = n zz >^(n) 

Conditioned on |A| = n, the distribution of A is a probability measure on Y„ 

which does not depend on £: 

M W (x) (z)x(z')x (dim A) 2 A (= Y 

M ^' [X) - zz'{zz' + l)...{zz' + n-l) ~a^~' A G Y " 

(recall that Y„ is the set of diagrams with n boxes). This means that the z 
measure M ZjZ / ; ^ is the mixture of the probability measures Mfl, with varying index 
n E Z + by means of the negative binomial distribution n zz ',£, see [B02], [B04]. 
For applications to harmonic analysis on the infinite symmetric group one needs 
the measures M y z z , and their scaling limits as n — > oo, but it turns out that the 
"mixed" measures M z>z i£ have much better properties, and the large n limit can be 
replaced, to a certain extent, by the £ — > 1 limit transition. This was the starting 
point of our paper [B02]. In the present paper we are dealing with the "mixed" 
measures only. 
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2. A BASIS IN THE £ 2 SPACE ON THE LATTICE AND THE MEIXNER POLYNOMIALS 

In this section we examine a nice orthonormal basis in the £ 2 space on the 1— 
dimensional lattice. The elements of this basis are eigenfunctions of a second order 
difference operator. They can be obtained from the classical Meixner polynomials 
via analytic continuation with respect to parameters. 

Throughout the section we will assume (unless otherwise stated) that parame- 
ters z, z' are in the principal series or in the complementary series but not in the 
degenerate series. In particular, z, z' are not integers. 

Consider the lattice of (proper) half-integers 

rpl _ j , 1 r _5_3_1135 1 
' 2 V ■ • J 2' 2' 2' 2' 2' 2>'"'J' 

Elements of Z' will be denoted by letters x, y. 

We introduce a family of functions on Tl depending on a parameter a € Z' and 
also on our basic parameters z, z', £: 



\T(z - a + ^)t(z' - a+ ±) J 
F(-z + a+\, -z' + a+\;x + a+l; 7 ^r) 



T(x + a+l) 



x e Z', (2.1) 



where £?; C; it;) is the Gauss hypergeometric function. 

Let us explain why this expression makes sense. Since, by convention, parameters 
z, z' do not take integral values, T(x + z + |) and T(x + z' + ^) have no singularities 
for x £ Z'. Moreover, the assumptions on (z, z') imply that 

r(a; + z + ±)T(x + z' + i) > 0, T(z - a + ±)r(z' - a + §) > 0, 

so that we can take the positive value of the square root in (2.1). Next, since £ € 
(0, 1), we have £/ (£ — 1) < 0, and as is well known, the function w — > F(A, B; C\ w) 
is well defined on the negative semi-axis w < 0. Finally, although F(A, B; C, w) is 
not defined at C = 0, -1, -2, . . . , the ratio F(A, 5; C, w)/T(C) is well defined for 
all C e C. 

Note also that the functions Va(£; z, z', £) are real-valued. Their origin will be 
explained below. 

Further, we introduce a second order difference operator D(z, z', £) on the lattice 
Z', depending on parameters z,z',£ and acting on functions f(x) (where x ranges 
over Z') as follows 

D(z,z',0f(x) = y/t(z + x+±)(z' + x + ±)f(x + l) 



+ 



^{ z + x -\){z> + x-\)f{x -l)-(x + £(z + z' + x)) f{x). 



Note that D(z,z',£) is a symmetric operator in £ 2 (Z'). 
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Proposition 2.1. The functions ip a (x; z, z' where a ranges over U , are eigen- 
functions of the operator D(z, z' ,(;), 

D(z, z 1 , £)1> a (x; z, z', = a(l - OlM*; z, z' , £). (2.2) 
Proof. This equation can be verified using the relation 

w(C - A){C - B)F(A, B- C + 1; w) - (1 - w)C{C - l)F{A, B;C-1; w) 
+C[C — 1 — (2C — A — B — l)w]F(A, B; C; w) = 

for the Gauss hypergeometric function, see, e.g., [Er, 2.8 (45)]. □ 

The next lemma provides us a convenient integral representation for the functions 

Lemma 2.2. For any A, B e C, M e Z, and £ e (0, 1) we have 
F(A, B; M + 1; ^r) _ r(-^ + l)r M/2 (l - Q B 



r(M + 1) T(-A + M + 1) 



7 s . (2.3) 



M 

Here £ G (0, 1) and {w} is an arbitrary simple contour which goes around the points 
and v 7 ? in the positive direction leaving l/v 7 ? outside. 

Comments. 1. The branch of the function (1 — y/^uj)^ 1 is specified by the con- 
vention that the argument of 1 — \^lo equals for real negative values of ui, and 

the same convention is used for the function ^1 — 

2. Like the Euler integral formula, formula (2.3) docs not make evident the 
symmetry A <-> B. 

3. The right-hand side of formula (2.3) makes sense for A = 1,2,..., when 
r(— A + 1) has a singularity. Then the whole expression can be understood, e.g., 
as the limit value as A approaches one of the points 1,2, ... . 

Proof. Since both sides of (2.3) are real-analytic functions of £ we may assume that 
£ is small enough. Then we may apply the binomial formula which gives 

r M/2 a - v^" 1 (i - ^) B "~ M 

_ ^2 ^2 + 1 ) fc ( j? )' ^(k+l-M)/2 Jk-l-M 



kill 

After integration only the terms with k — I + M survive. It follows that the right- 
hand side of (2.3) is equal to 

(l-C) 5 V T(-A + M + l + l) (B)i , 

T(-A + M + l) ^ T(l + M + l)ll ? ' 

V ; i>max(0,-M) V ; 
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We may replace the inequality I > max(0, —M) simply by I > because for negative 
integral values of M (when we have to start summation from I = —M), the terms 
with 1 = 0,..., —M — 1 automatically vanish due to the factor T(l + M + 1) in the 
denominator. Consequently, our expression is equal to 

(1 - Q B F(-A + 1 + M, B; M + 1; Q _ F(A, B; M + 1; ^) 
r(M + 1) ~~ T(M + 1) ' 

where we used the transformation formula [Er, 2.9 (4)]. □ 
Proposition 2.3. We have the following integral representations 

i>a(x;z,z',(,) 

'T(x + z+\)T{x + z' + |)V r(*'-aH 



r(z-a+i)r(z'-a+i) 



T(z> + x+\) 



(i-O- 



— 2: + a — 7 



M 



and 



V> (z; z, z', £)ip a (y; z, z', £) = ip ZtZ , (x, y) 



(2tt«) 2 



/ / (l-V^wi)' 



-2 +a- 



7 / 1 _ VI 



(wi) {uj 2 } 



(i - v^ 2 ) 



-z-\-a- 



W2 



(2.5) 



where 



Jv(x + z + 1)T^ + z> + i)r( y + z + i)r( y + z' + ±) 
(x ' y) = r(, + z- + i)r ( , + z + i) (2 - 6) 



_ffere eac/i contour is an arbitrary simple loop, oriented in positive direction, sur- 
rounding the points and and leaving l/v 7 ? outside. We also use the convention 
about the choice of argument as in Comment 1 to Lemma 2.2. 

Proof. Indeed, (2.4) immediately follows from (2.1) and (2.3). To prove (2.5) we 
multiply out the integral representation (2.4) for the first function and the same 
representation for the second function, but with z and z' interchanged. The trans- 
position z <-> z' in (2.4) is justified by the fact the initial formula (2.1) is symmetric 
with respect to z <-> z'. As a result of this trick the gamma prefactors involving a 
are completely cancelled out, and we obtain (2.5) □ 

Proposition 2.4. The functions ip a = tp a {x; z, z' , £), where a ranges over Tl , form 
an orthonormal basis in the Hilbert space l 2 (Z'). 
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Proof. From (2.4) it is not difficult to see that the function ip a (x; z, z', £) has expo- 
nential decay as x — > ±00. Indeed, depending on whether x goes to +00 or —00 we 
arrange the contour in such a way that \u\ > 1 or \u\ < 1, respectively. 

In particular, ip a (x;z,z',£) is square integrable. Since tp a is an eigenfunction of 
a symmetric difference operator whose coefficients have linear growth at ±00, and 
since to different indices a correspond different eigenvalues, we conclude that these 
functions are pairwise orthogonal in l 2 (U). 

Let us show that HV'all 2 = 1- Take (2.5) with x = y. Then the whole expression 
simplifies because (2.6) turns into 1. Next, in the double contour integral, we 
replace the variable u>2 by its inverse. We obtain 



1 1 



(Mx;z,z',()f = ^jj (l - V§*)~' + °~ 5 (l - VS"Z')'~°~* 



2 / u>i \ auii au>2 



To evaluate the squared norm we have to sum this expression over x € U . We split 
the sum into two parts according to the splitting l! = Il_ U We take as the 
contours concentric circles such that |wi| < \ in the sum over Z'_, and > \ui 2 \ 
in the sum over U + . This gives us 



^(^ a (x;z,z',e)) 5 

xei 1 



U>2 — U)\ U>i U>2 J J U>1 — Ul 2 (jJ\ U>2 

|wi|<|u)2| |u>i|>|u) 2 | 



with 

■ £ / ,— \ —z'+a 



x (1 - ^i^y^' 2 (1 - vt^y ~ a ~~ 2 1 



, ^2 ' 

Recall that both contours go in positive direction. 

Let us transform the second double-contour integral: keeping the second con- 
tour fixed we move the first contour inside the second contour. Then we obtain 
a double- contour integral which cancels the first double-contour integral, plus a 
single-contour integral arising from the residue of the function ui\ — ► (uji — u>2)~ 1 '- 

F(co,lu)— = — — <b jzr— — = 1. 



2m J v ' 'u 2 2m J (1 - y/lw)(w - VI) 

Thus, we have shown that the functions ip a form an orthonormal family in ^ 2 (Z'), 
and it remains to prove that this family is complete. For x € Z', let 5 X stand for the 
delta function at x. Since the functions 8 X form an orthonormal basis, it suffices to 
check that 



) 2 



But this follows from the previous claim and the symmetry a <-> x established in 
the next proposition. □ 
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Proposition 2.5. The following symmetry relation holds 
ip a (x; z, z', = ip x (a; -z, -z', £)■ 

Proof. Using the classical formula 

r(A + i)r(A-|) = 



COs(7rA) 

and the fact that both x + \ and a + \ arc integers we check that 

T{z + x+\)T{z' + x+\) _ T{-z + a+\)T(-z' + a + \) 
T{z-a+ \)T{z> -a+\)~ T{-z x + \)T{-z> -x+\)' 

Applying this to (2.1) and using another classical formula, 

F(A, B; C; w) = (1 - wf- A - B F(C - A,C - B;C; w), 

see [Er, (2.9.2)], we get the required relation. 

Another way to prove the proposition is to make a change of the variable in 
integral (2.4): 

, w - V£ 

LO I ► LO = — — . 

s/Zu-1 

This is an involutive transformation such that <-> an d 00 ^ Vvf- As is 
readily verified, it leads to transformation (a, x, z, z') — > (.t, a, — z, — z'). □ 

Corollary 2.6. TTie functions ip a — tjj a (x;z,z',^) satisfy the following three-term 
relation 

(1 - O^a = Y^Cz-a+i)(z'-a+i)Va-l 



- a - 5)(«' - a - 2) V'a+i + (-a + £{z + z' - a)) W 

Proof. Under symmetry x <-> a (Proposition 2.4), this turns into the formula stated 
in Proposition 2.1. Of course, a direct verification is also possible. □ 

The formulas of Proposition 2.1 and Corollary 2.6 show that the functions 
ip a (x; z, z' , £) possess the bispectrality property in the sense of [Gr]. 

Proposition 2.7. One more symmetry relation holds: 

Mx- z, z', - (-i) x+a V- (-*; -z, -z', 0, x,»e z'. 

Proof. This follows from the relation 

F(A,B;C;w) _ i_ c r(^-C+l)r(5-C + l) 

f(C) ~ w r(A)r(B) 

^ F(A-C+l,i?-C+l;2-C;H 

X f(2 = C) ' CG ' 



MEIXNER POLYNOMIALS AND RANDOM PARTITIONS 



13 



see [Er, 2.8 (19)]. Another way is to make a change of the variable, ui i— > 1/uj, in 
integral (2.4). □ 

In the remaining part of the section we will explain how the functions ip a are 
related to the Meixner polynomials. 

Let Z + = {0, 1,2,...}. To denote points of Z + we will use now the symbols x, y, 
because the letters x, y were already employed to denote points of U . Recall that 
the Meixner polynomials are the orthogonal polynomials with respect to the weight 
function 

M/ Mcixncr/~N _ (P)x£ X _ F(/3 + x)£ X , . 

W ™ {X) ~ x\ T((3)x\ ' Xe£+ ' (2J) 

on Z + , where (3 > and, as before, £ € (0, 1). Our notation for these polynomials is 
Wl n (x; [3, £). We use the same normalization of the polynomials as in the handbook 
[KS] (note that in [KS], our parameter £ is denoted as c). 
Set 

M n (x; P, = (-1)" ^ M ™(x), i£Z + , (2.8) 

where 

oo 

II2M • ; ft Oil 2 = E ^ & C)^ M r nCT (i). 

The factor (—1)™ is introduced for convenience: it will compensate the same factor 
in formula (2.10) below. 

Proposition 2.8. Drop the assumption that (z, z') is not in the degenerate series, 
and assume, just on the contrary, that z — N and z' — N+/3—1, where N = 1, 2, . . . 
and (3 > 0. Then expression (2.1) for the functions tp a { x : z , z ' ■, s ^ makes sense 
provided that the numbers 

x:=x + N-\, n:=N-a-\ (2.9) 

are in Z +; and in this notation we have 

ip a (x;z,z',£) = 9Jt„(:r;/3,£). 

Proof. We start with the expression of the Meixner polynomials through the Gauss 
hypergeometric function (see [KS, (1.9.1)]): 

m n (x;/3,0=F(-n,-x;p;^). 

Applying the transformation 

^ n ^ P ' W > T(p + n) r(l-6-n) ' " ' ' ' ' ' | 

we obtain 

_ (-l)"r(.x + l)r(/3) Fj-r^-p-n +hx + l-n:^ T ) 

T(f3 + n) \ £ J T(x + l-n) ■ i ■ ) 
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Although the first expression for the polynomials looks simpler than the second 
one, it turns out that only the second expression is suitable for our purposes. Note 
that (see [KS, (1.9.2)]) 

llJMn[ ' p ' m ~ r(/3)i> + i) • 

From the last two formulas and the definition of 2JL we obtain 



^ , l T{x + l)T{x + P )_ {i _ n)/2 _ 



(/3+2n)/2 



X 



F(-n, -(3 - n + l;x + 1 - n; ^r) 
T(x + l- n) 



Comparing this with (2.1) and taking into account (2.9) we get the required equal- 
ity. □ 

Thus, our functions ip a can be obtained from the Meixner polynomials by the 
following procedure: 

• We replace the initial polynomials Wl n by the functions Wl n ■ This step is quite 
clear: as a result we get functions which form an orthonormal basis in the £ 2 space 
on Z + with respect to the weight function 1. 

• Next, we make a change of the argument. Namely, we introduce an additional 
parameter N = 1, 2, . . . and we set x = x — N+^. Then we get orthogonal functions 
on the subset 

{-N+l-N+l-N+l,...}cZ', 

which exhausts the whole U in the limit as N goes to infinity. 

• Then we also need a change of the index. Namely, instead of n we have to take 
a = N — n — i. We cannot give a conceptual explanation of this transformation, it 
is dictated by the formulas. Again, the range of the possible values for a becomes 
larger together with N, and in the limit as N — > +oo we get the whole lattice U . 

• Finally, we make a (formal) analytic continuation in parameters N and f3, using 
an appropriate analytic expression for the Meixner polynomials (namely, (2.10)). 

We hope that this detailed explanation will help the reader to perceive the ana- 
lytic continuation arguments in section 3. 

Of course, instead of the lattice Z' we could equally well deal with the lattice 
Z, and then numerous "|" would disappear. However, dealing with the lattice Z' 
makes main formulas more symmetric. 

Remark 2.9. Note that the difference equation of Proposition 2.1 can be obtained 
via the procedure described above from the classical difference equation satisfied 
by the Meixner polynomials. This is precisely the way how we have obtained the 
difference operator D. Likewise, the three-term relation of Corollary 2.6 precisely 
corresponds to the classical three-term relation for the Meixner polynomials. 

3. The discrete hypergeometric kernel 

Let X be a countable set. By a point configuration in X we mean any subset 
X C X. Let Conf(X) be the set of all point configurations; this is a compact 
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space. Assume we are given a probability measure on Conf(X) so that we can 
speak about the random point configuration in X. The nth correlation function of 
our probability measure (where n — 1,2,...) is defined by 

p n (xi, ■ ■ ■ , x n ) = Probjthe random configuration contains x\, . . . , x n }, 

where x\, . . . ,x n are pairwise distinct points in X. The collection of all correlation 
functions determines the initial probability measure uniquely. 

We say that our probability measure is determinantal if there exists a function 
K (x, y) on X x X such that 

p n (xi, . . . ,x n ) = det [K(xi,Xj)]™ j=1 , n = l,2, ... (3.1) 

It is worth noting that if such a function K(x, y) exists, then it is not unique. 
Indeed, any "gauge transformation" of the form 

K(x,y)^^\K{x,y), (3.2) 

f{y) 

where / is a nonvanishing function on X, does not affect the determinants in the 
right-hand side of (3.1). 

Any function K(x,y) satisfying (3.1) will be called a correlation kernel of the 
initial determinantal measure. Two kernels giving the same system of correlation 
functions will be called equivalent. 

As in §2, we are dealing with the lattice U of (proper) half-integers. We split it 
into two parts, 1! = Z'_ U Z' + , where Z'_ consists of all negative half-integers and 
Z' + consists of all positive half-integers. For an arbitrary A € Y we set 

X(\) = {\ l -i + \ |i = l,2,...}cZ'. 

For instance, X_ (0) = U_. The set X (A) is sometimes called the Maya diagram of 
A, see, e.g. [M JD] . 

The correspondence A w 1(A) is a bijection between the Young diagrams A 
and those (infinite) subsets X C U for which the symmetric difference AA Z'_ is a 
finite set with equally many points in "L' + and Z'_ . Note that 

X(A') = -(Z'\X(A)). 

Using the correspondence A i— » X (A) we can interpret any probability measure 
M on Y as a probability measure on Conf(Z'). This makes it possible to speak 
about the correlation functions of M. Our goal is to compute them explicitly for 
the z-measures. 

Now we can state the main results of the paper. 

Theorem 3.1. Under the above correspondence between Young diagrams and Maya 
diagrams, any z-measure determines a determinantal measure on Conf(Z'). 
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Theorem 3.2. The correlation kernel of any z-measure M ZtZ /^ from the principal 
or complementary series can be written in the form 

K z , z ,^(x,y) = ^ ip a (x;z,z',£)ip a (y;z,z',t), x,yeZ', (3.3) 

where the functions ip a are defined in (2.1). 

Note that the series in the right-hand side is absolutely convergent. Indeed, 
since {ip a } is an orthonormal basis in £ 2 {Z') (Proposition 2.4), this follows from the 
fact that the series can be written as 

( S x,i>a)(lpa,Sy), 

where 5 X stands for the delta-function at point x on the lattice Z', and (■, ■) 
denotes the inner product in l 2 (Zl). 

Formula (3.3) simply means that K_ z , z ',^(x, y) is the matrix of the orthogonal 
projection operator in £ 2 (Z') whose range is the subspace spanned by the basis 
vectors i\) a with index a G Il + C Z' . 

Theorem 3.3. The correlation kernel (3.3) can also be written in the form 

Kz,z*,d x >y) = ( Pz,z'(x,y)Kz, z ',i(x,y) (3.4) 

where, as in (2.6), 



/T(x + z + l)T(x + z' + \)T{y + z + \)T{y + z' + \) 
^ AX ' V) = " r{x + z , + m y + z + l) (3 ' 5) 



and 
Kz,z',i{x,y) 



(2iri) 2 J J — 1 

{ui} {w} 



1 1 



x uj 1 2 uj 2 2 duJidu} 2 (3.6) 

where {oji} and {w 2 } are arbitrary simple contours satisfying the following three 
conditions: 

• both contours go around in positive direction; 

• the point £*/ 2 is in the interior of each of the contours while the point £ -1 / 2 
lies outside the contours; 

• the contour {wf 1 } is contained in the interior of the contour {w 2 } (equivalently, 
{uj^ 1 } * s contained in the interior of {ui\}). 

The kernels K_ ZtZ '.^(x,y) and K_ Z)Zt ^(x,y) are equivalent. Namely, they are re- 
lated by a "gauge transformation", 

Kz,z',dx,y) = { Z ' Z '\ X -, Kz,z',d x >y)> x >y e z '> 
jz,z' (y) 
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where 

/„(»)- . riX+ f + i) , (3-7) 

TTie kernel K_ ZiZ i^(x,y) can serve as a correlation kernel for the degenerate series 
as well. 

Proof of Theorems 3.1-3.3. We prove these three theorems simultaneously. Let 
Pn' Z ■ ■ ■ , x n) denote the n~point correlation function of M ZjZ >^. The proof 

splits into two parts. 

In the first part, we compute p"n' ,z '® for special values of the parameters corre- 
sponding to the degenerate series: z = N = 1,2,... and z' = N + (3 — 1, where 
(3 > 0. Here we use Proposition 1.4. We show that the formula 

pl z ' z ''t\ Xl , ...,x n ) = det||C x j)\i,j=i 

is valid (in particular, the values of the kernel in the right-hand size are well defined) 
when z — N, z' = z + [3—l, provided that N is so large that the numbers Xi + N— \ 
are nonnegative. Then we check that in that formula, the kernel K_ ZlZ >£ can be 
replaced by the kernel K_z,z'£- 

p^' z '^(x u ...,x n ) = det[K z . z ,. i (x l , Xj )]i tj=1 

In the second part, we extend the latter formula to other admissible values of 
parameters (z,z'). To do this we show that both sides are analytic functions in 
parameters (z,z',£). Moreover, these functions are of such a kind that they are 
uniquely determined by their values at points (z = N, z' = N + (3 — 1, £). 

We proceed to the detailed proof. 

Lemma 3.4. Let z = N = 1,2,... and z' = z + (3 — 1 with f3 > 0. Assume that 
x\, . . . ,x n lie in the subset Z + — N + \ C U , so that the points Xi := Xi + N — ^ 
are in Z+. 
Then 

p^'^ixu ...,x n ) = det [K^^ixi, i,-)] " J=1 , 

where 

JV-l 

K%»$* r (x, y)=J2 P> P> 0. Z+, 

m=0 

and the functions 9Jl m (x;/3, £) are defined in (2.8). 

Proof. According to Proposition 1.3, K^ ^ el (x,y) is the correlation kernel of the 
TV-point Meixner orthogonal polynomial ensemble with parameters (3 and £. 

On the other hand, let, as above, Y(N) denote the set of Young diagrams A with 
£(X) < N. Recall the bijective correspondence (1.5) 

X^X(X) = {x 1 ,...,x N } = {X 1 +N-l, Aa + JV-2,..., X N } 

between diagrams A € Y(N) and .ZV-point configurations in Z + . Comparing the 
definition of the infinite configuration X_ (A) C U with that of the A^-point config- 
uration X(X) we see that 

X(X) = (X{X) + N- |)nz+. 
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By Proposition 1.4, under this correspondence, the degenerate z-measure with 
parameters z — N, z' = N + (3 — 1, and £ turns to the TV-point Meixner ensemble 
with parameters (3 and £. This implies our claim. □ 

We take (3.3) as the definition of the kernel K_ z ,z',i{x,y). 

Lemma 3.5. Let z = N = 1,2, . . . and z' = z + (i — 1 with (3 > 0. ^4ssitme t/iat a; 
and y lie in the subset Z + — N + | c Z', so ffcai x := £ + AT — | and y := y + N — ^ 
are in Z+. 

T/ien expression (3.3) /or i/ie kernel K_ z . z '.^(x,y) is well defined and we have 
K z , z ,,s(x,y) = K™^(x,y). 

Proof. We have to prove that 

jv-i ^ 

£ ^ a (x; z, z', OMv; z, z', = S & $ & 

a6Z^ m=0 

We recall that the functions ip a (x; z, z', £) were defined under the assumption that 
both z, z' are not integers. However, as it can be seen from (2.1), each summand 
in the left-hand side of (3.8) makes sense under the hypotheses of the lemma. 

Set 

a(m) = N - m - §, m = 0, 1, . . . , N - 1. 
By Proposition 2.8, 

1pa(m)(x;Z,z',Z) = yR m {x;(3,C), Va(m) (V, Z, z' , = 9^m(j/; /?, 0) 

which implies that 

^ Va(a;; z, 0V>a(y; *, = £ & $ 0- ( 3 - 9 ) 



m=0 



Finally, observe that 
1 



T{z-a+\) 



a=AT+7T, N+ 



I w+ l r(iV-a+i) 



= 0. 



a=JV+i,JV+|, 



We conclude that the infinite sum in the left-hand side of (3.8) actually coincides 
with the finite sum in (3.9). □ 

Together with Lemma 3.4 this implies 

Corollary 3.6. Let z = N = 1, 2, . . . and z' = z + [3 — 1 with (i > 0. Assume that 
x\, . . . ,x„ lie in the subset Z + - TV + \ C II , so that the points Xt := x t + N - \ 
are in Z + . 
Then 

p ( n z - z '' 6) {x 1 , . . .,x n ) = det \K ZtZ , t t(xi, Xj)]" j=1 . 
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Lemma 3.7. Assume that 

• either (z, z') is not in the degenerate series and x,y G Tl are arbitrary 

• or z = N = 1, 2, ... , z' > N — 1, and both x, y are in Z + — TV + \ . 

Then the kernel K_ ZjZ / ^(x , y) of Theorem 3.3 is related to the kernel K_ z ,z',c(x,y) 
by equality (3.4). Equivalently, the kernels are related by the "gauge transformation" 
(3.2), 

where f z . z > is defined in (3.7). 

Proof. Let us start with expression (3.3) of the kernel K r r , f and let us replace 
each summand by its integral representation (2.5). It is convenient to set a — i = k 
so that as a ranges over U + , k ranges over {0, 1,2,...}. Then we obtain 

Kz,z',d x >v) = <Pz,z>(x,y) 

/ /('-^-•'(^yVv^-if " 



' (1 — v / € £j - , i)(l — V? w 2) \ k dwi du 2 



We can choose the contours {ooi} and {0J2} so that they are contained in the domain 
\u>\ > 1. Since the fractional-linear transformation 

preserves the unit circle |w| = 1 and maps its exterior |w| > 1 into its interior 
|w| < 1, we have on the product of the contours a bound of the form 



(1 - v^wi)(l - y/lw*) 



(wi - y/l){w 2 - V?) 



< q < 1. 



Therefore, we can interchange summation and integration and then sum the arising 
geometric progression in the integrand: 



£ 

fe=0 



1 ^ (l ^\ 

1 1 ) ^1^2 



(wi - ve)(w 2 - v 7 ?) ; " " (i-om-i) 



Then we obtain equality (3.4) with integral (3.6), as desired. Finally, we can re- 
lax the assumption on the contours: it suffices to assume that {w^ 1 } is strictly 
contained inside {^2}, as in the formulation of Theorem 3.3. 

It remains to show that (3.4) is equivalent to (3.10) According to (3.5) consider 
the expression 

1 T {x + z' + \)T{y + z+\) 



Vz.z'(x,y) 



Jr(x + z + i)r(x + z' + i)r( y + z + i)I\ y + z' + \) 
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Let us show that 

<Pz, Z '{x,y) fz, Z '(y) 

Indeed, l/<p z , z i has the form 

a(x)b{y) 
V ' a(x)b(x)a(y)b(y) ' 

and our hypotheses imply that a(x)b(x) and a(y)b(y) are real and strictly positive. 
We also have 

a(x) 



fz,z'(x) 



y/a(x)b(x) 
Therefore, we get 



fz.z'(x) _ a(x)y / a(y)b(y) _ a(x)a(y)b(y) 



fz,z'(y) y/a(x)b(x) a(y) y / a(x)b(x)a(y)b(y) a(y) 

= a ( x )b{y) = 1 

y/a(x)b(x)a(y)b(y) <Pz,z'{x,y)' 

□ 

Corollary 3.8. Let z = N = 1, 2, . . . and z' > N - 1. Then 

p} l * 1 * > «(ai,...,a; n ) = det ^^,5^, a:,-) (3.11) 

provided that all the points x\, . . . , x n G U lie in the subset Z + — N + \ CZ'. 

Proof. Indeed, this follows from Lemma 3.7 and Corollary 3.6. □ 

This completes the first part of the proof. Now we proceed to the second part. 

Lemma 3.9. (i) Fix an arbitrary set of Young diagrams V C Y. For any fixed 
admissible pair of parameters (z,z r ), the function 

a~ ^m^(a), 

\ET> 

which is initially defined on the interval (0, 1), can be extended to a holomorphic 
function in the unit disk |£| < 1. 

(ii) Consider the Taylor expansion of this function at £ = 0, 

00 

AeX> k=0 

Then the coefficients Gk,v{z, z') are polynomial functions in z, z' . That is, they are 
restrictions of polynomial functions to the set of admissible values (z,z'). 
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Proof, (i) Set V n = V n Y„. By the definition of M z . z ,.^ 

E M ^«( A ) = E ( E M Sw J 

AeX> n=o \Aex>„ / 

=(i-o"'f;f e m S( a )) 1££ ^- 

n=0 \.AeX>„ / 

Each interior sum is nonnegative and does not exceed 1. On the other hand, 



E frwfrOI = li - €1"' E {ZZ 'h < °°. IfK 1 - 

n=0 n=0 

This proves the first claim, 
(ii) By (1.11), 



E M x , z ,,t(x) = (i - o zz ' E E Wa^Oa^ 

AGX> n=0A€P„ 



/ dim A N 2 



V n! 



It follows that 



n=0 v ; AeX>„ V 

Since each V n is a finite set, this expression is a polynomial in z, z' . □ 

Now we can complete the proof of the theorems. Fix n and an arbitrary n-point 
subset X = {xi, . . . , x n } C Z', and regard p„' z '^(xi, . . . , x n ) as a function of 
parameters z, z',£. We want to show that equality (3.11) holds for any admissible 
(z, z'). Apply Lemma 3.11 to the set V of those diagrams A for which X_ (A) contains 
X, and observe that 



p^'^{x u ...,x n ) = Y,M z , z .,^)- 



Aex> 

It follows that pn' z '^(xi, . . . , x n ) is a real-analytic function of £ e (0, 1) which 
admits a holomorphic extension to the open unit disk |£| < 1. Moreover, the 
Taylor coefficients of this function depend on z, z' polynomially. 

On the other hand, from the expression (3.6) for the kernel K_ z ,z',^(x, y) it follows 
that this kernel (and hence the right-hand side of (3.11)) has the same property, 
with £ replaced by v 7 ?- 

Thus, both sides of (3.11) can be viewed as (restrictions of) holomorphic func- 
tions in t/£ with polynomial Taylor coefficients. Since the set 

{(z, z') | z is a large natural number N and z' > N — 1} 

is a set of uniqueness for polynomials in two variables, we conclude that equality 
(3.11) is true for any admissible (z,z'). 

This proves Theorem 3.1 and Theorem 3.3. Now, Theorem 3.2 follows from 
Theorem 3.3 and Lemma 3.7. □ 
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Proposition 3.10. Formula (3.3) for the kernel can also be written as 

y/I7^ i>-i(x)i>i(y)-i/>i(x)i>_i(y) 



x-y 



(3.12) 



Comment. The indeterminacy 0/0 arising on the diagonal x = y is resolved as 
follows. Observe that the defining analytic expression (2.1) for ip a {x) makes sense 
for any complex x sufficiently close to the lattice Z', so that we may view ipi ( • ) and 
"0_i( • ) as analytic functions in a neighborhood of Z' C C. Since the numerator in 
(3.13) is an analytic function in (x,y) vanishing on the diagonal x = y, it can be 
divided by x — y. Thus, the value of (3.12) on the diagonal can be computed, say, 
using the analytic expression (2.1) and the l'Hospital rule. 

Proof. Assume first x ^ y. Then it suffices to prove that 



{x-y) i>a{x)il> a {y) = -j— W>_i(x)V>i(y) - Vi(z)V>- 1 (</))• 

Recall the three-term relation from Corollary 2.6, which we can write as 

xip a {x) = A(a, a + l)ip a+1 (x) + A(a, a)ip a (x) + A(a, a - l)ip a _i(x) 

with appropriate coefficients A( ■ , ■). Using this and the similar relation for ytp a (y), 
and taking into account the symmetry relation 

A(a,a± 1) = A(a± l,a) 

(which follows from the explicit expression in Corollary 2.6), we readily get, after 
obvious cancellations, 

{x-y) ^a(x)ip a (y) = A(±,-±)(ip_i(x)ipi{y) ~ ipi(x)ip_i(y)). 

Since 

we are done. Notice that the infinite sums involved in this computation are conver- 
gent, because, for fixed x and y, ip a (x) and tp a (y) decay exponentially as a — > +oo: 
Indeed, by virtue of Proposition 2.5 this fact reduces to that pointed out in the 
beginning of proof of Proposition 2.4. 

To handle the case x = y we use the same trick as in Lemma 3.9 and the 
subsequent argument: the Taylor expansion at with respect to variable r\ := 

Specifically, let us regard ip a (x) = ip a (x; z, z', £) as a function in x e Z', a e Z' + , 
and 7] — From the integral representation (2.4) it is clear that this function is 
well defined as an analytic function in n ranging in the open unit disc \n\ < 1 in 
C. The same argument as above shows that this function decays exponentially as 
a — > +oo, uniformly on compact subsets of the disc. It follows that the series 

Kz,z',e{x,y) = Y tl>a,{x;z,z',£)il) a (y;z,z',t) (3.13) 

aGZ' + 
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is analytic in the same disc \rj\ < 1, too. 

On the other hand, from (2.1) it follows that if x + a > then ip a (x; z, z', £) is of 
order 0(i] x+a ) about r\ = 0. 2 Therefore, expanding the kernel in the Taylor series 
at r\ = 0, 

oo 
n=0 

we see that only finitely many terms in the series (3.13) contribute to any fixed 
coefficient F n (x, y; z, z'). Looking again at (2.1) we see that each coefficient can be 
written as 



F n (x,y; z, z') 

y/ (r(z + x + i)r(z' + x + i)r(z + y + %jr& + y + \) 

= T(x+\)T{y + \) 



G n (x,y;z,z'), 



where G n (x,y; z, z') is a rational function in (x,y). 

It follows that once we know the kernel out of the diagonal x = y we can 
extend it to the diagonal uniquely, by an obvious extension of the rational functions 
G n (x, y; z, z'). Finally, viewing the right-hand side of (3.12) as an analytic function 
in three variables x, y, and r/, it is readily checked that the recipe of extension 
suggested in the comment to the statement of the Proposition is the correct one. □ 

Remark 3.11. 1. The correlation functions of the z-measures M ZjZ /^ were first 
computed in [B02] in a different form: in that paper we dealt with another 
embedding of partitions into the set of lattice point configurations. The kernel 
K_ ZlZ '^(x,y) with x,y > coincides with one of the "blocks" of the kernel consid- 
ered in [B02]. The relation between both kernels is discussed in detail in [B05]. 
The proofs in [B02] and [B05] are very different from the arguments of the present 
section. 

2. Two other derivations of the kernel K r r i Ax, y) are given in Okounkov's pa- 
pers [Ok2] and [Okl]. In both these papers, the correlation functions are expressed 
through the vacuum state expectations of certain operators in the infinite wedge 
Fock space. A (substantial) difference between the methods of [Ok2] and [Okl] 
consists in the concrete choice of operators. The general formalism of Schur mea- 
sures presented in [Okl] is complemented by explicit computations in [BOk, §4]. 
One more derivation of the kernel K ■ y c(x. y) was recently suggested in [BOS]. 

3. In general, kernels of the form 

P(x)Q(y) - Q{*)P{y) 

x-y 

are called integrable kernels, in accordance to the terminology of [IIKS], [Del], [B]. 
In our case P and Q are expressed through the Gauss hypergeometric function, this 
is why we called K_ z ,z',£(x,y) the discrete hypergeometric kernel. 

2 It is worth noting that this claim is no longer true for negative x + a, because then the 
hypergeometric function in the numerator of (2.1) has a singularity compensated by the gamma 
function in the denominator, and the order at £ = has to be evaluated in a more sophisticated 
way. 
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4. The derivation of (3.12) from (3.3) is quite similar to the standard derivation 
of the Christoffel-Darboux formula for an arbitrary system of orthogonal polyno- 
mials. Since, as explained in §2, the functions ip a are closely related to the Meixner 
polynomials, this similarity is not surprising. 

5. Once we know that the functions ip a form an orthonormal basis (Proposition 
2.4), the series expression (3.3) for the kernel K_ z ,z',^{x, y) immediately implies that 
it is a projection kernel. This fact was first proved in [B05, §5] in a different way. 

6. The series representation (3.3) is equivalent to formula (3.16) in [Ok2]. A dou- 
ble contour representation of various correlation kernels related to Schur measures 
appeared earlier in [BOk]. 

7. Thus, almost all the results obtained in this section were already known. 
What is really new in our paper is the approach to their derivation based on the 
relationship to the Meixner polynomials. In [B07] we apply the same approach to 
a more complex (dynamical) model. 

8. One more novelty of the present work is appearance of the difference operator 
D; its importance becomes especially clear in the study of the dynamical model, 
sec [B07]. 

4. Remarks on a relationship to Krawtchouk polynomials 

There exists one more possible choice of basic parameters z, z', and £ leading to 
a family of probability measures on Y: Namely, parameters z, z' should be nonzero 
integers of opposite sign, while £ should be a negative real number (thus, instead 
of assuming £ £ (0, 1) we now require £ < 0). 

Indeed, let z — N and z' = —N', where N and N' are positive integers, and let 
£ < 0. Then the weight M Z>Z /^(A), as defined in (1.3), vanishes unless < N and 
£(X') < N', that is, A must be contained in the rectangle NxN'. For such diagrams 
A, we have (z)\ > while the sign of (z')\ equals (— 1)' A ' (see (1.2)). Since, the sign 
of £l A l also equals (-1) |A| , we have (z')a£ |A| > 0. Therefore, M z>z >^(\) > 0. The 
sum of all the weights is still equal to 1, so that we obtain an additional family of 
probability measures on Y. Let us call it the second degenerate series of z-measures. 
Its existence was pointed out in [B05, Example 1.6]. 

Let L be a positive integer and p £ (0, 1). The Krawtchouk weight function with 
parameters (p, L) is defined on the finite set {0, 1, ... , L} by 



The orthogonal polynomials with this weight are called the Krawtchouk polynomi- 
als, see [KS, §1.10]. Let us denote them as A n {x;p, L), where n is the degree of the 
polynomial. 

The next claim is a counterpart of Proposition 1.4 and can be checked directly: 

Proposition 4.1. Under correspondence (1.5), the z-measure of the second degen- 
erate series with parameters (z = N, z' — —N', £), where N,N' £ {1,2, . . .} and 
£ < 0, turns into the N -point Krawtchouk orthogonal polynomial ensemble with 
parameters 




x = 0, . . . , L. 



V = 



£-1' 



L = N + N' -1. 



Note that our assumption £ < implies < p < 1 . 
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The Krawtchouk polynomials &n(x;p, L) are close relatives of the Meixner poly- 
nomials 9Jl n (x; /?,£): both families of polynomials can be defined by the same ana- 
lytic expression involving the Gauss hypergeometric function, only the ranges of the 
parameters are different. The correspondence between the families can be formally 
written as follows: 

&n(x;p,L) = Wt n (x;-L,^ I ), 
see the very end of §1.10 in [KS]. 

Claim 4.2. In all arguments of the present paper that rely on the Meixner polyno- 
mials and the Meixner ensembles one could equally well use the Krawtchouk poly- 
nomials and the corresponding ensembles. 

Notice that any z-measure Mjv.-w.j of the second degenerate series can be 
written as a mixture of certain probability measures M^'_ N , living on the finite 
sets Y„, < n < NN', cf. Remark 1.5. Only now the "mixing distribution" on 
n's is not the negative binomial distribution but the ordinary binomial distribution 
with weights ( N f)p n {\ - p) NN '- n , where p = 

The measures M^_ N , admit a nice interpretation: Let Djv.jv' denote the rect- 
angular Young diagram with N rows and N' columns. Each standard tableau T of 
shape □jv.jV' can be viewed as a sequence of growing Young diagrams 

T = (0 = A<°>, A« . . . , AW . . . , A<™'> = n N , N ,), 

where A (n) S Y„. 

Proposition 4.3. Let {T} be the set of all standard tableaux of shape □jv.jv' 
equipped with the uniform probability measure. The push-forward of this measure 
under the projection T > A^"-* coincides with M^_ N , . 

In this form, the measures M^_ N , appeared in [PR]. A slightly different (but 
essentially equivalent) interpretation can be found in [B09, §5]. 

Notice that the measures M^, mentioned in Remark 1.5 can be obtained from 

the measures M^_ N , by analytic continuation in the parameters z, z': this ap- 
proach is developed in [B03]. 

Finally, notice that the material of this section is also related to the model 
considered in [GTW]. 
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